Magnetocapacitance of a graphene monolayer 
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Abstract 

We present a theoretical study of magnetocapacitance in a graphene monolayer at finite tem- 
perature taking into account the effects of disorder. The density of states (DOS) and magnetoca- 
pacitance found for graphene are compared to those found in standard two dimensional electron 
gas (2DEG) systems. The magnetic oscillations in DOS and magnetocapacitance are found to be 
enhanced and much more robust with respect to temperature damping in monolayer graphene in 
comparison with a 2DEG. Furthermore, we find that there is a it phase shift between magnetic 
oscillations in the two systems which can be attributed to Dirac electrons in graphene acquiring a 
Berry's phase as they traverse a closed path in a magnetic field. 
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I. INTRODUCTION 

Monolayer graphene is a gapless semiconductor with conical touching of electron and hole 
bands. Since the quasiparticles in a graphene monolayer obey the massless Dirac equation 
they act as massless Dirac Fermions. This difference in the nature of the quasiparticles in 
graphene from conventional 2DEG has given rise to a host of new and unusual phenomena 
such as anomalous quantum Hall effect and a n Berry phase Q, [J [J. Besides the funda- 
mental interest in understanding the electronic properties of graphene there are also serious 



suggestions that it can serve as the building block for nanoelectronic devices J4]. 

Among the most important tools in studying the electronic properties of solid state sys- 
tems are capacitance measurements. Capacitance measurements can be effectively used to 
probe the thermodynamic density of states of an electron system. It has been success- 
fully used toprobe the density of states in a conventional 2DEG systems and references 
therein In the present work, we determine the magnetocapacitance of a graphene 

monolayer. It is possible to carryout experiments, along the lines of those carried out on 
conventional 2DEG systems 0, 01 realized in semiconductor heterostructures, carbon nan- 
otubes and graphene nanoribbons[8[], that allow verification of the model presented here. If 
confirmed experimentally, our theoretical predictions provide another method of determining 
the density of states of a graphene monolayer in a magnetic field. 

In section II, we present the formulation of the problem. Section III contains the calcula- 
tion of the temperature dependent density of states and magnetocapacitance of a graphene 
monolayer, including the comparison with 2DEG and discussion of numerical results. The 
conclusions are in section IV. 



II. FORMULATION 

We consider two-dimensional Dirac electrons in graphene moving in the x-y-plane. The 
magnetic field (B) is applied along the z-direction perpendicular to the graphene plane. 
We employ the Landau gauge and write the vector potential as A = (0,Bx,0). The 
two-dimensional Dirac-like Hamiltonian l|, 2], [J for an electron in the Landau gauge is 
Hq = vj^a* .{— iHV + eA). Here a = {a x ,a y } are the Pauli matrices and vn charac- 
terizes the velocity of Dirac electron. The Landau level energy eigenvalues are given by 
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e n = hujgy/n, where u g = VD-y ^fp is the cyclotron frequency of Dirac electrons and n is 
an integer. The Landau level spectrum for Dirac electrons is significantly different from 
the spectrum for electrons in a conventional 2DEG which is given as e(n) = ftw c (n + 1/2), 
where uj c is the cyclotron frequency. This difference between the energy dispersion of Dirac 
electrons and conventional electrons is also reflected in the Density of States (DOS). The 
oscillations in the capacitance, in the presence of a magnetic field, are directly related 
to the density of states at the Fermi energy. In the presence of an external magnetic 
field, without taking into account any disorder effects, the DOS is represented as a se- 
ries of delta functions. In all practical situations, there is always some degree of disorder 
present in the system. The DOS delta functions broaden due to the presence of scat- 
tering centers and the DOS can be represented as a Gaussian, or a Lorentazian function 
to take this broadening into account [3] . Here, we first consider the Gaussian broaden- 



ing of DOS denoted by D G {e) given as D G (e) = ^y^exp 



(e-£ n ) 2 



2r 2 



, where T 



is the Gaussian distribution broadening width of zero shift. If we use the Poisson sum- 



: ~F(0) + ^TF(k) = / F(n)dn + 2^(-l) fc / F(n) cos[2Tikn}dn, the 

fc=l k=l 

summation over n in the DOS expression can be carried out by changing the variable 
as x = £ ~^Z" and applying the limit s ~ ep 3> T, where Bp = Vpfikp denotes the 
Fermi energy with Fermi wavenumber fcp = \f2im s and n s is the number density of the 
system. The first integral yields Dq = which is the DOS in the absence of mag- 



netic field and the second integral becomes proportional to / dx exp[— x 2 ] cos 
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ten as D G (e) = D I 1 - 2^(-l) fc exp 



The final result for the Gaussian broadened DOS can be writ- 
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k=l 



cos 



2nke 2 



Furthermore, follow- 



ing the same approach we can also determine Lorentazian broadened DOS D^e) 



t% 7 -i \i i ri with the result 



J D L (e) = J D Jl-2^(-l) fc exp 



k=l 



2-KkeT' 
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(1) 



We see that both forms of broadening yield the same numerical results 1(| in the regime 
s ~ Ef ^ r. For T ^> hujg, it is sufficient to retain only the first order term (k = 1 term) 



3 



since the contribution from the higher order terms is highly damped. Furthermore, the 
above result for the DOS is the same as obtained in Ref. 11 in the limit when energy gap 
A = [11]. 



III. MAGNETOCAPACITANCE OF A GRAPHENE MONOLAYER 



We consider a top-gated graphene device in which the capacitor is formed between the 
top gate and the graphene sheet. In order to deposit charge on the capacitor formed by the 
graphene sheet and the gate, the voltage source has to do both electrostatic and chemical 
work on the system. The effective capacitance is the geometric capacitance in series with 
"chemical capacitance". In an external magnetic field, the effective capacitance C(B) can 
be determined from the following relation [a, |6| 



1 



1 



1 



Dn 



D T (B) 



- 1 



(2) 



C(B) C e 2 D 

where Co = 1168 * 10~ 6 F m~ 2 is the capacitance at zero magnetic field and Dt{B) is 
the temperature dependent DOS at finite magnetic field that is determined from Dt{B) = 

o 

distribution function. In the above expression for Dt{B), we introduce a change of variable 
P(e — Ep) = s and apply the low temperature limit such that j3 <C Ep where (3 = -^p. As a 
result, it can be expressed as 



D T {B) = Do I 1 + exp 



2-KEpY 

H 2 Ug J 



cos 



2tte 2 f 

H 2 Ug J 



cos 



ds- 



cosh 2 [s/2] 



(3) 



The integration can be performed by using the following identity 12]: j dx -^^^ 

o 

-, with the result that the temperature dependent DOS is 



2f3 2 sinh(a7r/2/3) ' 



D T {B) = D \ 1 + 2- 



a 



sinh a 



cxp 
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where a = is a dimensionless parameter which expresses the characteristic temperature 
for damping of oscillations in Dt{B). 
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Following the same approach as given above for graphene monolayer, we can obtain the 



temperature dependent DOS for conventional 2DEG 



0, H, Q, E 



DAB) =DAl + 2^^- exp [-^ I cos 
sinh a, nw r 



9| with the result 
2tce p 



(5) 



where D s = ^ and a s = 

Now, we are in a position to compare the results for the temperature dependent DOS 
for monolayer, given by Eq.(4), with the result for 2DEG, shown in Eq. (5). The following 
differences in temperature dependent DOS in the two systems are found: 

i) -The argument of the oscillatory cosine term, COs[^f - 7T] responsible for the 7r-Berry 
phase shift, in 2DEG depends on the ratio of the Fermi energy to the cyclotron frequency 
while for monolayer graphene the argument of the cosine depends on the square of the ratio 
of the Fermi energy to the corresponding cyclotron frequency without an extra n. 

ii) -The exponential term, exp[— responsible for the exponential decay of the ampli- 
tude of the oscillatory DOS for 2DEG depends on the ratio of T to the cyclotron frequency 
u c . From Eq.(4), in monolayer graphene the exponential term depends on the ratio of the 
product of T and Ef to the square of the corresponding cyclotron frequency u g . The magni- 
tude of the argument of the exponential term is 1.84 in GaAs 2DEG while 0.22 in Graphene 
at fixed value of the magnetic field (1 Tesla) and disorder (0.5 meV). Therefore, we ex- 
pect that for experimentally relevant parameters as considered here, disorder will cause less 
damping of DOS oscillations in a graphene monolayer compared to a standard 2DEG. 

iii) - the temperature dependence of the oscillation amplitude of DOS is given by the 
factor lim . °, , = 2ae~ a . The amplitude of the oscillatory DOS decays exponentially as 

a ,_> 0O sinn^a) 

the temperature is raised. From this temperature dependent coefficient the characteristic 
temperature for damping of oscillations in DOS, To, can be identified such that ksTn = 
l^j^j with kp being the Fermi wavenumber, equal to \J2ixn s , while in a standard 2DEG[s|, [f], 
the corresponding expression is ksT s = = . Here m* = 0.068 m e in GaAs and 
m e is the usual electron mass. These characteristic temperatures determine the robustness 
with respect to temperature of magnetic oscillations in DOS. The relative magnitudes of 
the characteristic temperatures for damping in 2DEG and graphene monolayer is ^ = 
the ratio of the characteristic temperatures is equal to the ratio of the corresponding Fermi 
velocities in the two systems, where v p = is the Fermi velocity in a 2DEG. Typically, the 



Z, 
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Fermi velocity of 2DEG systems (~ 10 5 m/s) is an order of magnitude smaller than that of 
Graphene (~ 10 6 m/s). An estimate of this ratio of temperatures is given by ^- = = 
- 0.24 for GaAs, For Ge is ~ 0.03 and ~ 0.015 for Si with a fixed value of the 



number density ( n s = 3.2 x 10 15 m™ 2 ). For comparison for an equal temperature scale the 
effective mass should be ~ 0.016 m e , which is smaller than that for most 2DEG systems 
such as Si, Ge, GaAs. 

From the discussion of the DOS expressions in the two systems given above, we expect 
that the oscillations in the magnetocapacitance will be affected less by temperature in a 
monolayer graphene compared to a 2DEG (GaAs is considered for comparison in the nu- 
merical results) due to the higher damping temperature in the monolayer. On varying the 
magnetic field, changes in the DOS are reflected by oscillations of the magnetocapacitance. 
Since magnetocapacitance is directly related to the DOS, the exponential decay in the am- 
plitude of the magnetocapacitance oscillations with temperature and disorder is expected. 
Before we present results for the magnetocapacitance, the dimensionless DOS (Dt(B) / D ) 
for the two systems is shown in Fig.(l). The DOS is plotted as function of the inverse mag- 
netic field. The dimensionless magnetic field Bq/B is introduced through the characteristic 
magnetic field B = hk 2 F . B = 13.25 T for n s = 3.2 x 10 15 m~ 2 . The temperature is 1.2 
K and the disorder energy T = 0.5 meV in Fig.(l). The magnetic oscillations, with 7r-Berry 
phase difference, are clearly visible in the DOS of the two systems . Furthermore, the os- 
cillation in the DOS in the graphene monolayer are found to be enhanced and less damped 
with temperature compared to those in a 2DEG. In order to study the effects of temperature 
on the magnetocapacitance in the monolayer and 2DEG, we show in Figs. (2, 3) the tem- 
perature dependent magnetocapacitance as a function of the magnetic field at two different 
temperatures: 1.5 K and 22 K. In these figures, T is held constant at 0.5 meV. We observe, as 
expected from the above discussion, that oscillations in magnetocapacitance are less damped 
by temperature and there is a 7r-Berry phase difference in graphene monolayer compared 
with the same oscillations found in a 2DEG. For the parameters 5|, |6|, [7j considered here and 
with constant disorder, the oscillations found in the magnetocapacitance in a 2DEG become 
completely damped at 22 K whereas they are found to persist upto a temperature of 60 K 
in a graphene monolayer. Therefore, we have found that magnetocapacitance oscillations 
are more robust with respect to temperature damping in a graphene monolayer compared 
to a 2DEG. Furthermore, these magnetic oscillations have a 7r-Berry phase difference in the 
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two systems. 



IV. CONCLUSIONS 

We have investigated the finite temperature density of states and magnetocapacitance in a 
graphene monolayer taking into account finite temperature as well as disorder. The results 
obtained are compared to those found for conventional 2DEG systems. We have found 
that magnetic oscillations in the DOS and magnetocapacitance in a graphene monolayer 
are enhanced, are more robust against temperature and have a 7r-Berry phase difference 
compared to a 2DEG. 
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